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One-dimensional spin- 2 XY models are among the few many-body prob- 
lems for which a lot of exact results can be obtained. This is due to the 
fact that the Hamiltonian of the system can be rewritten via Jordan- Wigner 
transformation in terms of non- interacting fermions [jl|. In spite of this there 
are rather limited results for time- dependent correlations between longitu- 
dinal spin components. More difficult to deal with are the random versions 
of the model. Apart from academic importance such results are used in the 
description of a large number of quasi-one-dimensional ferroelectric and mag- 
netic compounds such as C s{Hi_xDx)2P0i, PhHi_x D^-PO^, PrCl^ etc. In 
our previous paper we suggested a new numerical approach for evalua- 
tion of thermodynamical properties and both equal-time and different-time 
spin correlation functions for spin-| non- uniform anisotropic XY chain in 
transverse field with the Hamiltonian 

N N-l 

H=T. + E (^r^^j+i + (1) 

This method allows to study the random versions of the model (1) assuming 
that JJ^, J™, ... vary randomly from site to site according 

to the given probability distribution. The purpose of the present paper is to 
describe in more detail the proposed approach and to apply it for examining 
the infiuence of disorder on transverse dynamical susceptibility of spin-| 
Ising chain in random transverse field. Such results, up to the best authors' 
knowledge, are obtained for the first time. 

In order to obtain the equilibrium statistical properties of the system in 
question (1) one should first perform Jordan- Wigner transformation from the 
raising and lowering operators = sj±isj to Fermi operators Ci = Cj = 
n{=i(-2<)s7,j = 2,...,iV, c+ = s+,c+ = m-\(-2<)4,j = 2,...,iV and 



then to diagonalize the obtained quadratic form by a canonical transforma- 
tion rjk = T,f=i{gkjCj + hkjC^) with the outcome H = T,k=i^k{VkVk - i)- 
The elementary excitation spectrum Aj. and the coefficients g^n = {^kn + 
^kn)/'^: hkn = {^kn " *fcn)/2 are determined from 

N N 

JXX I jyy JXX I jyy 

A — o X _i_ i ~^ i X _i_ i-i X 

JXX _ jyy JXX _ jyy 
Bij = — — ^ ' ^ ' (2) 

Eqs. (2) reduce to standard problems 

N 

A^$,„ = ^$,,[(A-B)(A + B)],„, 

TV 

A^^,„ = ^ *,,[(A + B)(A - B)],„ (3) 

i=i 

foT N X N five diagonal banded matrices (A =f B) ( A ± B) . 

Next it is necessary to express the quantities of interest through A^, 
^kj, ^fcj- One immediatelly notes that the evaluation of partition function 
Z{/3, N) = llk=i 2 cosh ^ involves all A|. Thus, the knowledge of eigenval- 
ues of matrices (A=pB)(A±B) permits one to obtain the thermodynamical 
properties of the model in question (1). Somewhat more complicated is the 
investigation of spin correlation functions. Reformulating spin operators in 
terms of operators cpf = Cj' ± cj (they are the known linear combinations of 
operators r]k, r]^, namely, (p+ = Em=i ^mjiVm + Vm), = Em=i "^mjiv^ - 
rjj), that is, = m-\(99+^-)(^+/2, = Ui~=\{vtVn)vJ /m, s| = 
—ifj'ipj/2, noting that the calculation of spin correlation function reduces to 
exploiting Wick-Bloch-de Dominicis theorem and calculating the elementary 



contractions 



2 



I / N I x-^ cosh (i Apt — 



p=i cosh ^ 

< >= E'^P.^P™ cosh ^ 

^ sinh(-zApt + ^) 



p=i ^™ cosh^ 

^ coshfiA t - ^) 

< >= - E J^M, ' , (4) 

p=i cosh 

one finds that ^fem for all A^ are involved in final expressions for 

both equal-time and different-time spin correlation functions. Therefore, the 
study of equilibrium statistical properties may be performed numerically in 
a following way: first, to solve the eigenvalues and eigenvectors problem for 
matrix (A + B)(A — B) (3) obtaining in result A^ and '^kji and second, 
having '^^j and A^ = -y/A|, to find from Eq.(2); thermodynamics and 
spin correlation functions are expressed via the sought quantities. It should 
be stressed that here appears only the eigenvalues and eigenvectors problem 
for N X N five diagonal banded matrix that is the remarkable peculiarity 
of the system in question (1). This is a key difference of our approach in 
comparison with other finite-chain calculations (reported e.g. in 0) where 
the numerical diagonalization of the Hamiltonian (1) in the Hilbert space of 
dimension 2^ is performed (see also 0, |^). 

We end up with presenting the results of numerical calculation of trans- 
verse dynamical susceptibility Xzz{^,^) = Y.n^^^"' lo^ dte^^'^~^^''^^h < [s|(t), 
Sj+n] > fo^ spin- 1 Ising chain ( JJ^ = J, Jj^ = 0) in random transverse field 
rij with probability distribution density p(...,fij, ...) = Ylf=i[xS{ilj) + (1 — 
x)6{Qj — Q)]. Acting as described above one finds that 4 < s|(t)s|^„ >=< 



Vtit)'^j it)vt+nVj+n > = < '^t^j >< ^i+n^j+n > " < ¥?+(t)</?++„ >< 
it)Vj+n > + < Vtit)y^j+n >< Vjit)vt+n >, < > = < s|+„(-t)s| >. 

Computing then numerically the elementary contractions (4) involved in the 
spin correlation functions < sf^Qit) s^qq^^ >, < sf^Q^^^—t) s^qq > for the ran- 
dom chain of 200 spins with J = —2, Q = 1 ai j3 = 20, performing nu- 
merically the integration over t and summation over n, and averaging over 
the realization (typically few hundred) we have obtained the results some of 
which are presented in Fig.l. They agree with generally known results for 
perfect (non-random) cases of Ising model when x = 1 and Ising model in 
transverse field when x = derived for the first time with periodic boundary 
conditions imposed in As the concentration of the sites with transverse 
field Q increases Xzz{ce,uj) rebuilds from the Ising-like type behaviour to 
the behaviour inherent to Ising model in transverse field. The divergence 
of ReXzz{ce,u!) at u = 2 becomes finite, decreases, loses its symmetry and 
develops into a smooth curve in the limiting = 0. The (5-singularity 

at = 2 of lmXzz{ce, uj) transforms into few sharp cusps, loses its symmetry 
moving slightly to the right and also develops into a smooth curve at x = 0. 
The study of the dependence of presented results on number of spins in the 
chain, on upper limit in integration over t, on value of e, on number of terms 
in sum over n and on number of random realizations were performed care- 
fully. It should be noted that the obtained susceptibilities exhibits a lot of 
structure (especially at small x) that requires further study. 

Since a complete solution for some pair time-dependent correlation func- 
tions even in perfect spin-| XY chains has never been found, the numerical 
calculations seem to be of much use in deriving new results here. Some rig- 
orous results on longitudinal time-dependent correlation function were de- 



rived numerically for perfect one- dimensional spin-i transverse Ising model 
at critical field and isotropic XY model in zero field at /5 = oo in 0, how- 
ever, within quite different approach. The suggested method, that gives 
exact results for finite spin-| XY chains, appears to be of great use in study 
of dynamical properties and random versions of systems considered. The 
derived conclusions may be useful for interpretation of observable data in 
dynamical experiments on quasi-one- dimensional compounds that can be 
described in frames of model (1). 

The present paper was presented as a poster at International Conference 
on Magnetism 1994 (22-26 August, 1994, Warsaw, Poland) (EP.227) g]. 
The authors would like to express their gratitude to Prof. L.L.Gongalves 
and other participants of ICM'94 for interest in this paper and discussions. 
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Fig. 1. Frequency-dependent real and imaginary parts of transverse sus- 
ceptibility Xzz{^,^) for Ising chain in random transverse field at ae = for 
different values of concentration x; e — 0.01. 



Fig. 1. Frequency-dependent real and imaginary parts of transverse sus- 
ceptibility Xzz{^,^) for Ising chain in random transverse field at ae = for 
different values of concentration x; e — 0.01. 
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